Extension the Noether's theorem to Lagrangian formulation with 

nonlocality 

Zaixing Huang 

State Key Laboratory of Mechanics and Control of Mechanical Structures 

Nanjing University of Aeronautics and Astronautics 

Yudao Street 29, Nanjing, 210016, P R China 

E-mail: huangzx@nuaa.edu.cn 



^D 



43 



> 



February 25, 2013 



Abstract 



A Lagrangian formulation with nonlocality is investigated in this paper. The nonlocality of the Lagrangian 
Oh| is introduced by a new nonlocal argument that is defined as a nonlocal residual satisfying the zero mean condi- 

tion. The nonlocal Euler-Lagrangian equation is derived from the Hamilton's principle. The Noether's theorem 
is extended to this Lagrangian formulation with nonlocality. With the help of the extended Noether's theorem, 
the conservation laws relevant to energy, linear momentum, angular momentum and the Eshelby tensor are de- 
termined in the nonlocal elasticity associated with the mechanically based constitutive model. The results show 
that the conservation laws exist only in the form of the integral over the whole domain occupied by body. The 
localization of the conservation laws is discussed in detail. We demonstrate that not every conservation law cor- 
responds to a local equilibrium equation. Only when the nonlocal residual of conservation current exists, can a 
conservation law be transformed into a local equilibrium equation by localization. 

Key words: nonlocal Euler-Lagrange equation, Noether's theorem, conservation law, nonlocal elasticity, me- 
chanically based constitutive model 



££ ■ 1 Introduction 

(N 

Nonlocal elasticity has developed into an important branch of continuum mechanics since the first pioneer studies 
in the last 60-70s [1, 2, 3, 4, 5]. This theory and the extension of it have been applied to various topics in engineer- 
ing [5, 6, 7, 8]. So far, constitutive models in the nonlocal elasticity can be categorized into three sorts : Eringen's 
constitutive model [4, 7], peridynamic constitutive model [9, 10, 11] and mechanically based constitutive model 
[12, 13]. The postulation of the Eringen' constitutive model consists in that the stress at a point not only depends 
on strain of the point but also on strains of all points within body. So stress is equal to the convolution integral of 
strain over the domain occupied by the body. Due to this fact, the mixed boundary value problem is ill-posed in 
the nonlocal elasticity associated with the Eringen's constitutive model [14], except the attenuating kernel being 
equipped some ad hoc futures, for example, it is assumed to be the Green function of the differential operator V 2 . 
Under the latter assumption, the integral-type constitutive model reduces to the so-called implicit gradient model. 

In the peridynamic constitutive model [9, 10, 11], ones introduce the internal long-range body force to represent 
the interactions within body, but forsake the conception of stress and strain. An integral operator is used to for- 
mulate the constitutive relation between the internal long-range body force and displacements. Since there are no 
assumptions made on the differentiability of displacements in the motion equation of the peridynamic nonlocal 
elasticity (peridynamics), this theory is suitable to study phenomena with discontinuities and fragmentation. In the 
peridynamics, no boundary conditions appear as there are no spatial derivatives. In order to solve the boundary 
value problems, some unconventional boundary conditions have been prescribed [11, 15]. Unlike the conventional 
boundary conditions, they are imposed on a boundary layer with non-zero volumetric measure, as opposed to a 
geometric boundary in the strict mathematical meanings. This inconsistency may give rise to some difficulties in 
the problems with complicated boundary conditions [12, 13]. 

The mechanically based constitutive model (MBCM) can be regarded as a fusion between the Eringen's consti- 
tutive model and the peridynamic constitutive model [16]. This model retains the conception of stress and strain 
but introduces meanwhile the internal long-range body force characterizing the interactions between non-adjacent 



particles. The internal long-range body force linearly or nonlinearly depends on the relative displacements be- 
tween interacting particles within body; while the stress-strain relation is still characterized by the conventional 
constitutive equation, e.g., the Hooke's law. Stress and the internal long-range body force are independent on each 
other. If a nonlocal elasticity is concerned with MBCM, we refer it as to the nonlocal elasticity associated with 
MBCM. In this theory, all types of boundary conditions are the same as that in the classical elasticity. From the 
view of this point, it is of advantage to adopt MBCM in nonloacal elasticity. 

Paola et al [12, 13] firstly advanced the linear theory of the nonlocal elasticity associated with MBCM. They called 
it the mechanically based model of nonlocal elasticity. In this theory, the internal long-range body force at a par- 
ticle is linearly dependent of the relative displacements between the particle and other particles within body. If 
taking a suitable nonlocal kernel, this internal long-range body force can be also obtained by linearization of the 
the peridynamic constitutive model. 

Recently, the Lagrangian formulation has been proposed for the mechanically based model of nonlocal elasticity 
[17]. The relevant energy-momentum tensor is given under a simplest case in which the nonlocal kernel degrades 
into a constant equal to the reciprocal of the volume of body. The energy-momentum tensor is of fundamental im- 
portance because it represents the configurational force on a defect (e.g., vacancy, inclusion, dislocation and crack 
etc) in solids; while in the absence of defects, it embodies a conservation law [18]. In order to determine, under a 
more general case, the energy-momentum tensor and the relevant conservation laws, it is necessary to extend the 
Noether's theorem to the nonlocal elasticity associated with MBCM. 

The study on the nonlocal form of the Noether's theorem can be traced back to Edelen [19]. He [19, 20, 21] inves- 
tigated the reformulation of the Noether theorem in the nonlocal theory based on a general theoretical framework. 
In this framework, the nonlocal argument specified by Edelen [19] is too general in form to take account of the 
constraint of physical laws. Later, Edelen [22] simplified the nonlocal argument as a linear integral operator on 
the field variable. Under this case, Vukobrat and Kuzmannovic [23]) addressed the conservation laws in the non- 
local elasticity associated with the Eringen's constitutive model. Recently, Lazar and Kirchner [24] discussed the 
energy-momentum tensor and relevant configurational forces in the nonlocal theory. They issued some interesting 
results on the interaction between dislocation and disclination. 

At present, few results are known on the form of the Noether's theorem and relevant conservation laws in the non- 
local elasticity associated with MBCM. So the objective of current work is to clarify these subjects. The paper is 
organized as follows: we start in Section 2 by using a nonlinear integral operator to define the nonlocal argument. 
Based on this nonlocal argument, we develop a Lagrange formulation with nonlocality. The Noether's theorem is 
extended into the new Lagrangian formulation in Section 3. According to the extended the Noether's theorem, in 
Section 4 we investigate the conservation laws in the nonlocal elasticity associated MBCM in which the internal 
long-range body force is nonlinearly dependent of the relative displacements between particles. In section 5, the 
localization of the conservation laws is discussed. The relevant nonlocal residuals are determined. We close this 
paper in Section 6 by making some concluding remarks. 

Notation: A compact notation is used, with boldface letters being vectors or tensors. The index rules and sum- 
mation convention are adopted. Latin indices have range 1, 2, 3; while Greek indices run from to 3. Partial 
derivatives with respect to coordinates are represented as a j = da/dx k . Partial derivative with respect to time is 
denoted by an upper dot, i.e., a = da/dt. The norm | • | is defined as y/ {•){•). For example, \r\ — yjr-r. Other 
symbols will be introduced in the text where they appear for the first time. 

2 Lagrangian formulation with nonlocality 

A continuum occupies the domain Q. in the three-dimensional Euclidean space. Let every particle in the continuum 
be referred to by the orthogonal Cartesian coordinates x = {x',x 2 ,x 3 } specifying its position in £2, and let (p = 
(p(t,x) denote a field variable defined on £2. Depending on circumstances, <p is a scalar, vector or tensor. Let (h\q>) 
represent the nonlocal argument of (p. We define it as 

(%) = (p{t,x) I h{x,y, \cp(t,x) - (p(t, y)|)dv(y) - / h{x,y, \cp(t,x) - (p(t,y)\)(p(t,y)dv(y), (1) 

Jo. Jq. 



where h(x,y, \(p(t,x) — <p(f,y)|) is called the nonlocal kernel. Let r = (p(t,x) — (p(t,y). Eq. (1) can be abbreviated 
to 

(h\(p) = (p(t,x) I h{x,y,\r\)dv(y)- / h(x,y,\r\)(p(t,y)dv(y) 
Jo. Jn 

= I Kx,y, \r\)[(p{t,x) - (p(t,y)]dv(y) 
Ja 

h(x,y,\r\)rdv(y). (2) 



Different from the case recently reported by Huang [17], the nonlocal kernel depends not only on x and y but also 
on the field variables. Therefore, (h\(p) is a nonlinear integral operator with respect to (p. We enforce the nonlocal 
kernel to fulfil the symmetry: h(x,y, \r\) — h(y,x, \r\). Under this symmetry, it is easy to verify that (h\q>) satisfies 
the zero mean condition: 

r (h\(p)dv(x)=0. (3) 



la. 

Due to Eq.(3), the nonlocal argument determined by Eq. (1) is in essence different from the definition given by 
Edelen [19, 22]. The latter fails to satisfy the zero mean condition. 



LetL = L(t,x,(p,<p,(p t, (h\(p))(k = 1,2,3) denote the Lagrangian '. So the action functional of <p can be written 

as 

A[(p}= L(t,x,(p,f,(p Jc ,(h\(p))dv(x)dt. (4) 

Jt Ja 

Suppose <p(f ,x), h(x,y, \r\) and L are suitably smooth functions. In order to determine the variation SA[(p], we 
firstly prove the two identities. Using Eq. (2), we have 

6{h\q>) = f{[8h(x,y,\r\)]r + h(x,y,\r\)8r}dv(y)= f(h+\r\^-)8rdv(y) 
Jo. Jn a\r\ 

= 8cp(t,x) f (h + \rAdv(y)- f (h + \rA8cp(t,y)dv(y) 
Jn o\r\ Ja a\r\ 

= (g\8(p). (5) 

[ \l/(t,x)(h\<p)dv(x) = f f y/(f,x)/z(x,y,|r|)[<p(f,x)-<p(f,y)]dv(y)dv(x) 
J a JnJa 

(p(t,x)h(x,y, |r|)[v/(f,x) - \j/(t ,y)]dv(y)dv(x) 



nJo. 

<p(t,x)(h\\i/)dv(x). (6) 

a 

Eq. (6) is valid for any continuous function h, \j/ and (p. For example, we may use g in Eq. (5) to replace h, where 

g reads 

. dh , . .. . dh(x,y, I r I) 

g = h+ r =h(x,y, r \)+\r\ ' J ' u . (7) 

d\r\ d\r\ 

Clearly, g is symmetric with respect to x and y. By means of Eq. (5) and Eq.(6), we can calculate the variation of 

A[q>]: 

sm = ri ( ^ 5<p+ ^^ + ^^ + 4M 5</zi<p>)dv(x)df 

h f'i f dL 

dv(x)+ / / ^ — n k 8(pds(x)dt, (8) 

t Jt Jdad<p,k 



a d(p 



where dQ. is the boundary surface of Q. and n\ denotes the unit normal vector on d£l. In Eq. (8), a shortened form 
similar to Eq.(l) is used, 

' gfrj, I r |)dv(y) - / g(x,y, \ r |)_-L-dv(y). (9) 



d(h\<p)' d(h\<p)Jn 6y "" u w; Ja (,y,JU u d(h\<p) 



1 If necessary, the nonlocal arguments (h\(j>) and {h\ip^) may be conveniently inserted into L. But in this case, the boundary conditions will 
become too complicated to solve in mathematics. 



Let dQ. = dQ.\ U dQ.2, dQ.i n dQ.2 = 0. On dQi, <p takes a given value <p. Then, the boundary condition on dQ\ 
reads 

9\d^=9- ( 10 ) 

At the initial and terminal time, we have 

<Pl = 90, <Pl fl =9i- (11) 

Due to Eq.(lO), 8(p =0on dQi. Similarly, 8(p = at the initial and terminal time because of Eq.(ll). Thus, 
Eq.(8) reduces to 

8A[q>]= [" [[^-U^)-(4^)k + (g\J¥h)]$(pdv(x)dt+ f l [ ^n k 8(pAs{x)6d. (12) 
L J Jt Jndcp dt y dcp ; v <?<p/' XSf d(h\(p) n w J t0 Jdn 2 d(p, k w 

In terms of the Hamilton's principle, we have 8A [<p] = 0. So the fundamental lemma of variation yields the below 

results: 

Euler-Lagrange equation: 

d , dL . . dL . dL . . dL . 

d7 ( ^ )+ W ,ft "^ = ^ l ^%) ) - (3) 

Natural boundary condition: 

dL 



-n k 



= 0. (14) 

d& 2 



d<p 

Eq.(13) is also referred to as the nonlocal Euler-Lagrange equation. The right-side term of Eq.(13) is the nonlocal 
term, called the nonlocal traction. If h is independent of r, we have g~h, and then Eq.(13) will reduce to the case 
in [17]. 

Using the symmetry of g(x,y, \r\), we easily verify the equality below: 

Hm)L si "' y - ]r l)d,(,)dv(x) = Uj { ^ - n i»W) M > )M * ] - <15) 

by interchanging x and y. As thus, the integral of Eq.(9) over Q. leads to 

L {8l dm )Mx) =°> (16) 

which shows that the nonlocal traction is a nonlocal residual automatically satisfying the zero mean condition. Due 
to Eq.(16), the integral of Eq.(13) over Q. has the same expression as the integral of the Euler-Lagrangian equation 
not concerned with nonlocality. 

In physics, if Eq. (13) represents the equation of motion, the nonlocal traction may be interpreted as an internal 
long-range body force applied on a particle at x by other particles within body. Every particle is subjected to such 
a force. In terms of the action and reaction law, the sum of all such forces must be zero. Therefore, Eq.(16) is just 
an embodiment of the action and reaction law. 

3 Extension of the Noether's theorem to the Lagrangian formulation with 
nonlocality 

For convenience, time can be treated as a coordinate so as to form a four-dimensional position vector x with 
coordinates x = {x"} = {x° ,x l ,x 2 ,x 3 } defined on £2, where x° = t and Q. = [f ,? 1 ] Ufl. Thus, the Lagrangian 
function can be abbreviated as L = L(x,(p,(py, (<p)), (7 = 0, 1,2,3), and the nonlocal Euler-lagrange equation (see 
Eq.(13)) can be rewritten as 

. dL . dL . . dL . 

W ,r_ ^ = ^ l ^%) ) - (7) 

Consider infinitesimal transformations of group 

x 1 — > x' = x + c5x. (18) 

(p(x) ^ <p'(x') = (p(x)+A(p(x). (19) 

y ^ y' = y + 5y. (20) 

<p(y) — ► 9'(y') = <p(y)+A(p(y). (2i) 



Here, A<p is different from the variation 8(p. The latter is defined as 

c5<p = <p'(x)-<p(x). (22) 

Therefore, between A<p and 8(p, there exists the relation below: 

A<p = 8(p + <p. r <5x r . (23) 

The action functional A [<p] is said to be symmetry if it is form-invariant with respect to the infinitesimal transfor- 
mations (18) - (21), i.e., 

f L[x>,( P >(x>),( P > r (x>),(h\( P ')}dv(x>) = /.L[i,p(i),p, r (i),(%)]dv(i). (24) 

Jo! " J a 

In the left integral, x' now represents merely a dummy variable of integration and can therefore be relabeled x. But 
there remains a change in the domain of integration, so Eq. (24) becomes 

/ L[x,(p'(x),(p'(x),(h\(p')]dv(x) = [ L[x,<Kx),<p, r (x),(%)]dv(x). (25) 

It should be noted that (h\(p') in the left-hand term of Eq. (25) is an integral defined on D.', but the integral domain 
of (h\q>) in the right-hand term is £2. For (h\(p'), we have 

(%') = <p'(x) / h(x,y, |r'|)dv(y) - / h(x,y, \r'\)(p' (y)dv(y). (26) 

J a' J a' 

By the transport theorem [25], Eq. (26) leads to 

<%') = ^(x){^Mx,y,k1)dv(y) + ^ a[Mx '^ l)5/] dv(y)} 

d[h(x,yy\) (p'(y)8y* 
a in ay 



Afry, |r'|)<p'(y)dv(y) - /_ "VKWM v,», J dv(y) . (27) 



Expanding h(x,y, |r'|) at r to the first-order term yields 



= Mx,y,kl) + ^^^^(x)- 5 <p(y)]. (28) 



Inserting Eq. (22) and (28) in (27) and omitting the higher-order terms, we have 

<%') = <%) + {g \ 8(P ) + J a ^,yM)r8/ ] dv{yy (29) 

After Eq. (22) and (29)are substituted into L[x, (p'(x),(p'y(x), (h\(p')], it becomes 



L[±y(i),^(±),<%'>] = L[x, ( p + 8( P ,( P , r +(8( P ), r ,(h\(p) + (g\8(p) + j J^M0)I^l Av{y)] 

= L{x,<p,<p iV {h\<p)) + 8L, (30) 

where 

Cr dL s dL fS . dL .. . s . f d[h(x,y,\r\)r8y k ] J , ,, 

5L -^ + ^ 5 ^ + ww) M5<p)+ L w dv(y)] - (31) 

By Eq. (30), we have 

I L[x,(p'(x),(p'(x),(h\(p')}dv(x) = f {L[x,<Kx),<p, y (x),(%)]+SL}dv(x). (32) 

JQ! " JQ! 

Applying the transport theorem to the right-hand of Eq. (32) and omitting the higher-order terms, we have 

f L[x,<p'(x),<p' (x),(%')]dv(x) = / L[±,p(i),p, y (i),<%)]dv(i) + f 8Ldv(x) + f (L8x?) , y .dv(x). (33) 

JQ' " ' JQ JQ JQ 



Substituting Eq. (33) into (25) leads to 



L 



<5Ldv(x) + / (L8x r ) y .dv(x) = 0. (34) 



Using Eq. (31), we have 

/ Sidv(x) = / ' df / 8lAv(x) 
Jn Jt Q Ja 



'i f r dL e dL ,. , dL r . .. . /■ d[/t(x,y,|r|)r5;y* 



[ ^-W ,r+<g| a<* 

d r „ „ /■ <?L 



^-W' +w 3^ >1 * <,)+ jkW ,f) '' d * 



W 5 *7/^ y ' |r| >5<*M dvW1 < lv( * ) ' <35) 



By Eq. (17), Eq. (35) reduces to 



Substituting Eq. (36) into (34) leads to 

f A [^-S<p+LSxr-8xrJ a rh{x,y, |r|)^|^dv(y)], r dv(x) = 0. (37) 

If Eq. (18) and (19) belong to a finite Lie group, 8x? and A<p can be represented as [26] 

Sx? = e a xl A(p = e a <P a , (38) 

where e" is an infinitesimal parameter independent of the space-time coordinates. X% and <P a are the generators 
of Lie group. By using Eq. (38) and (23), Eq. (37) becomes 

L { -k*° + W-1**'- w mf) M " ))S ' - ;t^«>>*<*> = °- < 3! » 

Thus, the Noether theorem is extended to the nonlocal Lagrangian formulation with nonlocality. Let 
where J a is called the conservation current (or Noether current). Thus, Eq. (39) can be abbreviated as 



■111 



J aj dv(x)=0, or / 7> y da(x)=0, (41) 



where n y is an unit normal vector on d£l. Eq.(41) shows that total conservation current /„ on Ci. is conserved. 
However, Eq. (41) will cease to be valid if its integral domain Ci is replaced by a subdomain of A. This is because 
the nonlocal argument is defined on the whole domain £2, and it can not be altered in the deductive process from 
Eq.(24) to (41). In other words, if we use any v (v= [f ,? 1 ] Uv, v C £2) to replace Q. in Eq. (24) but keep the integral 
domain of the nonlocal argument fixed, it will be impossible to derive Eq. (41) from (24) due to needing to change 
the order of integrals through the interchange of x and y. Therefore, the conservation current is conserved on Q. as a 
whole, and yet when exactly the same statement is made for a subdomain of Q. it is no longer valid. As a result, we 
:,y = Ra (x) from the localization of Eq. (41), rather than Jay ■■ 



can merely derive /„ _ y = R a (x) from the localization of Eq. (41), rather than /« y = 0. R a (x) is called the nonlocal 



residual of the conservation current. It should satisfy the zero mean condition so that the integral of J a , y = Ra (x) 
over Ci can return to Eq. (41). Huang put forward a representation of the nonlocal residual automatically satisfying 
the zero mean condition [27]. 

For the convenience, in subsequent sections we call the Eq. (41) the conservation law, while J„ y = /?«(x) is 
referred to as the local equilibrium equation. 



4 Conservation laws in nonlocal elasticity 

Consider a linear, homogenous, nonlocal elastic body free of external body forces. For this body, the Lagrangian 
function is assumed to take the following form: 

1 1 1 

L = -pum - -QjkiUijUkj - -{h\Ui)uu (42) 

where m, denotes the elastic displacement field, p and Qm are the mass density and elastic tensor. The term 
(h\ui)Ui/2 represents the internal long-range action potential. 

Using m instead of <p in (13), and then inserting Eq. (42) in (13), we have 

ptii+(K\ut) = Qju Uk,ij , (43) 

where the K is represented as 

K = h(x,y,\r i \)+ 1 -\r l \ dh{x ^f ll) , n = Ul (x)- Ui (y). (44) 

2 d\rt\ 

It is easy to see that Eq. (43) characterizes the motion equation of a nonlocal elasticity associated with MBCM. 
This theory is an extension of the mechanically based model of nonlocal elasticity [12, 13, 17]. If h is supposed 
to be independent of |r,|, then K — h. Thus, Eq. (43) will reduce to the motion equation in the mechanicall based 
model of nonlocal elasticity, see [12, 13, 17]. 

In order to find the conservation laws corresponding to Eq. (43), we firstly calculate the conservation current 
according to Eq. (42). Let <p = u k . Then Eq. (40) becomes 

Jay=j. *ta + {[£- / r t h(x,y,\ri\) dv(y)]g y|U - 3 u k ^}X m - (45) 

o«t,y Ja o{h\ui) "^ OU k y 

It should be remembered that, in the convention of this paper, Latin indices take 1, 2 and 3; while Greek indices 
run from to 3. Since uq is null, <P k o = 0. Using Eq. (42) and (45), and noticing u kt o = u k = du k /dt, we have 

Jay,y = -j- [p«t*ta + (L+ - / h(x,y, \n\)r k u k (t ,y)dv(y))X 0a - pu k (u k X 0a + u k:S X sa )] 



[C k sijUij^> ktt -(L+- / h(x,y, \n\)r k u k (t ,y)dv(y))X sa - C ksij Uij(u k Xo a +u k}[ Xi a )] !S . (46) 



1 
2Jn 

By Hooke's law a ks = C kS jjUij, Eq. (46) is rewritten as 

J ay,y = -^ [ [puk'5>ka + {L+- h(x,y, \ri\)r k u k (t, y)dv(y))X 0a - pu k (u k X 0a + u ktS X sa )} 

- [<y ks $>ka-(L+- /!(x,y,|r,D^M A .(f,y)dv(y))X, a -a fa (M^ 0a + M M X /a )] j ,. (47) 

2 in 

In the classical elasticity, Fletcher proved the completeness of conservation laws under the infinitesimal transfor- 
mations below [28]: 

1 1 — > t' = t + e(vt + co) , xt i — ► x' t = xi + e(vxi + Eij k xjb k + ci), 

u t \ — > u • = u t + e (- vuj + e ijk Ujb k + e ijk xjc k +d t ), (48) 

where e is a infinitesimal parameter. V,ai,bi,c a and <f, are arbitrary real constants. In terms of Eq. (48), we turn 
now to investigating the concrete forms of the conservation laws under four typical transformations. 

4.1 Case 1: t' — t + e, x' k =Xk, u' k = iik 

The transformations above are equivalent to taking Xqq = l,X sa = Q^ ka = in Eq. (38). Under this case, Eq. 
(47) reduces to 

Substituting Eq. (49) into (41) yields 

-(L+- / h(x,y,\r i \)r k u k (t,y)dv{y)-pu k u k ) + (o ks u k ), s ]dv(x)=0. (50) 

n dt 2Jn 



In terms of Eq. (3), it is easy to verify that 

/ / h(x,y,\r i \)r k u k (t,y)dv(y)dv(x) = / u k (t,y) / h(x,y,\r { \)r k dv(x)dv(y) 
JaJa J a J a 

= [ u k (t,y) [ (h\u k )dv(x)dv(y) = 0. (51) 

Jn Jn 

Substituting Eq. (51) into (50) yields 

[—(L-pii k u k ) + (o ks u k ), s ]dv(x) = 0. (52) 

a df 



By Eq. (54) and the divergence theorem, Eq. (52) is rewritten as 

/• 1 1 1 f 

■ {^pu k u k +-Ci jk iUijU k } + -{h\u k )u k )dv(x)- I o ks u k n s da(x) = 0. (53) 

Jo. 2 z z j^a 



d 
d7. 



Eq. (52) corresponds to the conservation of energy that shows total energy on £2, including the sum of kinetic 
energy, elastic potential energy and the internal long-range action potential energy, is equal to work done by 
external traction. However, because the internal long-range interactions give rise to energy transferring among 
different parts of body, the same form as Eq. (52) or (53) is no longer valid for the subdomain of £2. -This is just a 
intrinsic character solely processed by the nonlocal theory. 

4.2 Case 2:t' = t, x' k = x k , u' k = u k + £k, £\ = £2 = £3 = £ 

The transformations above represent the rigid body translations. To satisfy these transformations, we set Xy a = 
0,<£>sjfc = 8 sk . Thus, Eq. (47) reduces to 

Jky,y = fa(P a k) - (Ofa)„s- (54) 

Substituting Eq. (54) into (41) leads to 

— / pu k dv(x)- / a ks n s dv{x) = 0, (55) 

dt Jo. J da 

which is the integral representation of the conservation law of linear momentum. Although Eq. (55) has the same 
form as the relevant conservation law in the classical elasticity, it can not be transformed into the same differential 
equation due to the nonlocality. 

4.3 Case 3: t' — t, x' k = x k , u' k = u k + £kijXi£j, £\ = £2 = £3 = £ 

The transformations above characterize the rigid body rotations. Under this case, we have X a n = 0,<t> k j = £ k tjXi. 
So Eq. (47) leads to 

Jky,y = -r;{p£kjiXiU}) - (£kjiXi<Jj S ), s - (56) 



Substituting Eq. (56) into (41) gives 

d 
dtJn 



/ pe kji XiUjdv(x) - / s k jiXi<Jj S n s dv(x) = 0. (57) 

JO. JdQ. 



Therefore, the conservation of total angular momentum on £2 is associated with the invariance of the action func- 
tional under the rigid body rotation. 

4.4 Case 4: t' — t, x' k = Xk + £ k , u' k = u^ £\ = £2 = £3 = £ 

The transformations above correspond to the coordinate translations that are identical to setting Xy = 8jj,Xo a — 
0,X a0 = 0,4> ta = 0. Thus, Eq. (47) reduces to 

hy,y = --T.iPUjU.ij;) + [(L+- h(x,y, \n\)r k u k {t,y)dv{y))8 jk + OijU lk \j- (58) 

Substituting Eq. (58) into (41) and using (51), we have 

— I puiU iyk dv{x) - (Ln k + G ij u Lk n j )dv(x)=0. (59) 

df Jn Jdn 



Let Tn — LSjic + Oijiiik. Tjk is the so-called Eshelby tensor. In the classical elasticity, the Eshelby tensor is 
independent of the integral path. However, the same conclusion is no longer available in nonlocal elasticity. This 
is because Eq. (59) holds only on Q. as a whole. For any v C £2, it ceases to be valid. Eq. (59) represents the 
conservation laws relevant to the Eshelby tensor. 

In the formal four-dimensional space of Cl, Eq. (53) and (59) can be combined in a unified expression with the 
help of the energy-momentum tensor. This expression has been given in [17]. 

5 Localization of conservation laws and nonlocal residuals 

In previous sections, we have pointed out that, if the nonlocality is concerned, the conservation laws are valid 
only on whole domain occupied by body but fail on the local domain of the body. Therefore, the localization 
of a conservation law is bound to result in the nonlocal residual of conservation current appearing in the local 
equilibrium equation, in which, the divergence of the conservation current is equal to the nonlocal residual rather 
than zero. Thus, by localization we can give the local equilibrium equations respectively corresponding to Eq. 
(52), (55), (57) and (59). They are formally represented as follows: 

— (L-pu k u k ) + (o ks u k ) :S = E. (60) 

-7-{pUlc)-(Oks),s=Pk- (61) 

— (pe kji XiUj) - (e k jiXiGj S )j = M k . (62) 

— {pinuik) - {L8 jk + CTijU Lk )j = J k . (63) 

Here, E denotes the nonlocal residual, P k the nonlocal residual of linear momentum, M k the nonlocal residual of 
angular momentum and J k the nonlocal residual of the Eshelby tensor. All nonlocal residuals should be satisfy the 
zero-mean codition, i.e., 

[ (E,P k ,M k ,f k )dv(x)=0, (64) 

Ja 

so that the integrals of Eq. (60), (61), (62) and (63) over Q. can return to Eq. (52), (55), (57) and (59), respectively. 

Now, we turn to determining the nonlocal residuals and the existence of Eq. (60) - (63). 

5.1 Nonlocal residual of energy 

The conservation laws are firmly associated with the characters of the LagrangianL(f,x,M,', m,-, u^ k , (h\ui)). In fact, 
we have noticed that if L does not depend explicitly upon f, Eq. (52) can be determined by calculating total 
derivation of L with respect to f in the integral sign of the action functional A [<p] . Thus, if directly evaluating the 
derivation of L with respect to f, we have 

dL dL dL dL dL dihlui) 

~a7 = 3~~ Ui + JT~ Ui + J> M '> + 3/; i \ a, 

dt out out au lk a{h\Ui) At 

r dL d dL dL , . d dL . , dL . , dL ... 

= b -- ^ ;)- " ;i — W M; + ^( ;)-"'') + (;)— u ^- k + lm^^- (65) 

out dt out oUi k at out duj k d\h\Ui) 

Applying Eq. (13) to (65) gives 

dL d . dL . , . dL . , dL . ... . . dL . . 

df dt dut du[ k d{h\ufi d{h\Ui) 



Eq. (66) can be rewritten as 



d dL . dL . , dL .... , dL 

-77 {L- -j-ui) - (-j Ui).k = -w,, A g\ui) - g 7577-1 , K (67) 

df dut dui k d\h\Ui) d\h\Ui) 



By using Eq. (42), Eq. (67) leads to 

_d 

df 
Compared Eq. (68) with (60), we have 



(L-pu t u i ) + (o ik u i ) k = Ui(t,x) / g(x,y,\ri\)ui{t,y)dv{x)-Ui(t,x) / g(x,y,\n\)ui(t,y)dv(y). (68) 

Ja Jo. 



E = iii(t,x) / g(x,y,\ri\)ui(t,y)dv(x)-Ui(t,x) / g(x,y,\n\)ui(t,y)dv(y). (69) 

Jo. Ja 

So far, E has been determined. Clearly, it satisfies Eq. (64). 



5.2 Nonlocal residuals of linear and angular momentum 

Eq. (61) characterizes the equilibrium of linear momentum. So it is identical to Eq. (43). Due to this fact, it is 
easy to see that/\ = — (K\u k ), i.e., 

P k = f K(x,y,\ ri \)u k (t,y)dv(y)-u k (t,x) f jr(x,y,|iv|)dv(x). (70) 

Jn Jn 

P k can be used to represent M k . In order to verify this point, we firstly rewritten Eq. (62) as follows: 

d 

M k = —(p£ k jiXiUj)-{£ k jiXiO js ) tS 

= £kjiXi[—{puj)-(a js ). s ]. (71) 

Inserting Eq. (61) in (71) gives 

M k = e kji XiPj. (72) 

In a general case, it is easy to see that e k pXiPj always violates the zero-mean condition. As a result, M k given by 
Eq. (72) does not satisfy Eq. (64). On the other hand, M k should follow Eq. (64) so as to ensure the integrals of 
Eq. (62) over Q. returning to Eq. (57). The paradox caused by Eq. (72) shows that it is impossible, under a general 
case, to transform Eq. (57) into the form of Eq. (62) by localization. Eq. (62) is inaccessible. 

However, if the nonlocal kernel takes the form of the central pair potential, £ k jtXiPj can be consistent with the 
zero-mean condition. Under this circumstance, (h\u k ) is simplified into (h\x k ), and P k can be represented as 

P k = [ Jc(x,y, \n\)y k dv{y) -x k f Jc(x,y, |r,|)dv(x). (73) 

Jn Jn 

Substituting Eq. (73) into (72) yields 



M k = £ kji xi / jc(x,y, |r,|)y/dv(y), (74) 

Jn 



which is equivalent to the expression below: 



M = xx / K(x,y,|r 1 -|)ydv(y). (75) 

Jn 



Since both x and y are the dummy variables of integration in the below, we have 

/ Mdv(x) = / / jc(x,y,|r,|)xxydv(y)dv(x) 
Jo, JnJn 

Jf(y,x, |r;|)yxxdv(x)dv(y) 



nJo. 



Jf(x,y, \n\)x x ydv(y)dv(x) = - / Mdv(x). (76) 

laJa J a 

Eq. (76) leads to 

/ Mdv(x) = 0, i.e., / M^dv(x) = 0. (77) 

Jn Jn 

So far, we have demonstrated that M k satisfies the zero-mean condition when the long-range body force is governed 
by a central pair potential. 

5.3 Nonlocal residual of the Eshelby tensor 

The conservation law relevant to the Eshelby tensor attributes to the invariance of the action functional under the 
coordinate translations. Therefore, both the Lagrangian L and the nonlocal kernel h are bound to be explicitly 
independent of x k . As thus, we have 

dL dL dL dL ., . . 

L k = 3— M(,fc+3^«(,fc+3 ui,jk+ -w, I x {h\Ui), k 

out out ouij o{n\U[) 

r dL d dL dL d .dL , dL , dL 

= b -- ^ ;)--;r- )j]"^ + ^-(j-«a) + (^— «.",* j + 57rrT *"')■*■ (78) 

out at aui aui.j df out ouij a{n\Ui) 



Applying Eq. (13) to (78) leads to 



d dL dL dL dL 

at dui diij.j ' a{h\Ui) d(h\ui) 



Substituting Eq. (42) into (79) yields 



Eq. (80) is rewritten as 



L,k = -r{puiU Lk ) - (aijU Lk )j + Ui(h\ui) :k - (g|«j) «,-,*. (80) 



— {puiUit) - {L8 kj + OijU Lk )j = (g\ui)u Lk - Ui{h\ui) )k . (81) 



Comparison between Eq. (80) and Eq. (63) will gives 

h = (g\ui)u Lk - Ui(h\ui), k . (82) 

Clearly, J k determined by Eq. (82) fails to agree with the zero-mean condition. - This means that the nonlocal 
residual of the Eshelby tensor does not exist. Therefore, Eq. (81) is not a local equilibrium equation corresponding 
to Eq. (59). In other word, Eq. (59) can not be transformed into a local form similar to Eq. (63). The conservation 
law relevant to the Eshelby tensor exists only in the form of integration. 

6 Conclusions 

In this paper, we extend the definition of the nonlocal argument by introducing a nonlinear nonlocal kernel de- 
pending not only on the spatial coordinates but also on the field variables. The extended definition retains the zero 
mean character of the nonlocal argument. It is this character to distinguish the nonlocal Lagrangian formulation 
developed in this paper from other nonlocal variational theories. 

On the basis of the extended nonlocal argument, a Lagrangian formulation with nonlocality is established. The 
nonlocal Euler-Lagrange equation is derived from the Hamilton's principle. Accompanied with this equation, the 
nonlocal traction appears in the form of the nonlocal residual satisfying the zero mean condition automatically. 
- This is an obvious difference between the new theory and the other theories. Physically, the nonlocal traction 
represents the long-range interactions within body. Therefore, the zero-mean character of the nonlocal traction is 
just a embodiment of the action and reaction law. 

The Noether's theorem is extended to the Lagrangian formulation with nonlocality. The result shows that conser- 
vation law exists only in the form of the integral over the whole domain occupied by body. The local equilibrium 
equation derived from the localization of the conservation law is equal to the nonlocal residual of the conservation 
current, provided it exists, rather than zero like the case in the variational theories without nonlocality. In physics, 
the presence of the nonlocal residual attributes to the conservation current transferring, caused by the internal long - 
range interactions, among different parts within body. 

A Lagrangian including the nonlocal argument is advanced in the quadratic form. The motion equation derived 
from this Lagrangian is consistent with the nonlocal elasticity associated with MBCM. In this theoretical framwork, 
we use the extended Noether's theorem to determine the conservation laws relevant to energy, linear momentum, 
angular momentum and the Eshelby tensor. The localization of these conservation laws are discussed in detail. 
We demonstrate that local equilibrium equation of energy and of linear momentum can be respectively derived 
from the corresponding conservation laws by localization, but the conservation law relevant to the Eshelby tensor 
can not be transformed into a local form by localization. So no local equilibrium equation relevant to the Eshelby 
tensor exists in the nonlocal elasticity associated with MBCM. The nonlocal residual of energy and nonlocal resid- 
ual of linear momentum have been determined, respectively. They are consistent in mathematical form with the 
representation of nonlocal residual given in [27]. 

In general, there is no local equilibrium equation corresponding to the conservation law of angular momentum. 
However, if the nonlocal kernel takes the form of the central pair potential, the local equilibrium equation of angu- 
lar momentum will occur in the nonlocal elasticity associated with MBCM. Under this case, the internal long-range 
interaction manifests itself as a central force field. 

We therefore conclude that, in the nonlocal elasticity associated with MBCM, not every conservation law cor- 
responds to a local equilibrium equation. Only when the nonlocal residual of conservation current exists, can a 



conservation law be transformed into a local equilibrium equation by localization. The results in this paper imply 
that he existence of the local equilibrium equation is to some degree influenced by the nonlocal kernel, However, 
that is a problem needing further exploration. 

Acknowledgements 

The support of the National Nature Science Foundation of China through the Grant No. GAA12012 is gratefully 
acknowledged. . 

References 

[I] Kroner, E. Elasticity theory of materials with long-range cohesive forces. International Journal of Solids and 

Structures, 13, 731-742 (1967). 

[2] Krumhansl, J. A. Some considerations of the relation between solid state physics and generalized continuum 
mechanics, in Mechanics of Generalized Continua, pp298-31 1, ed., E Kroner, Springer, New York, 1968. 

[3] Edelen, D.G.B., Green, A.E. and Laws, N. Nonlocal continuum mechanics. Archive for Rational Mechanics 
and Analysis, 43, 36-44 (1971). 

[4] Eringen, A.C. and Edelen, D.G.B. On nonlocal elasticity. International Journal of Engineerng Science, 10, 

233-248(1972). 

[5] Kunin, LA. Elastic Media with Microstructure (I and II), Springer, Berlin, 1982, 1983. 

[6] Rogula, D. Nonlocal Theory of Material Media, Springer, New York, 1982. 

[7] Eringen,A.C. Nonlocal Continuum Field Theories, Springer, New York, 2002. 

[8] Bazant, Z.P. and Jirasek, M. Nonlocal integral formulations of plasticity and damage: survey of progress. 
ASCE Journal Engineering Mechanics, 128, 1 1 19-1 149 (2002). 

[9] Silling, S.A., Reformulation of elasticity theory for discontinuities and long-rang forces. Journal of the Me- 
chanics and Physics of Solids, 48, 175-209 (2000). 

[10] Emmrich, E. and Weckner, O. Analysis and numerical approximation of an integro-differential equation 
modeling non-local effects in linear elasticity. Mathematics and Mechanics of Solids, 12, 363-384 (2007). 

[II] Silling, S.A. and Lehoucq, R.B. Peridynamic theory of solid mechanics. Advances in Applied mechanics, 44, 
73-168 (2010). 

[12] Di Paola, M., Failla, G. and Zingales, M. Physically-based approach to the mechanics of strong nonlocal 
elasticity theory. Journal of Elasticity, 97, 103-130 (2009). 

[13] Di Paola, M., Pirrotta, A. and Zingales, M. Mechanically-based approach to nonlocal elasticity theory: vari- 
ational principles. International Journal of Solids and Structures, 45, 539-548 (2010). 

[14] Atkinson, C. Crack problems in nonlocal elasticity. Archive of Mechanics, 32, 597-614 (1980). 

[15] Aksoylu, B. and Mengesha, T Results on nonlocal boundary value problems. Numerical Functional Analysis 
and Optimization, 31, 1301-1317 (2010). 

[16] Zingales, M., Di Paola, M. and Inzerillo G. The finite element method for the mechanically based model of 
non-local continuum. International Journal for Numerical Methods in Engineering, 86, 1558-1576 (201 1). 

[17] Zaixing Huang. Lagrangian formulations with nonlocal residual-based arguments. Annals of Solid and Struc- 
tural Mechanics, 2, 99-105 (201 1). 

[18] Maugin, GA. Material Inhomogeneities in elasticity, Chapman & Hall, Londen, 1993. 

[19] Edelen, D.G.B. Nonlocal Variations and Local Invariance of Fields, Elservier, New York, 1969. 

[20] Edelen, D.G.B. Invariance theory for nonlocal variations principle-II: weak invariance and global conserva- 
tion laws. International Journal of Engineering Science, 9, 801-814 (1971). 



[21] Edelen, D.G.B. Invariance theory for nonlocal variations principle-Ill: Strong and absolute invariance. Inter- 
national Journal of Engineering Science, 9, 815-829 (1971). 

[22] Edelen, D.G.B. Nonlocal Field Theories, in Continuum physics (IV), pp75-204, ed., A.C. Eringen, Academic 
Press, New York, 1976. 

[23] Vukobrat, M. and Kuzmanovic, D. Conservation laws in nonlocal elasticity. Acta Mechanics, 92, 1-8 (1992). 

[24] Lazar, M. and Kirchner, H. O. K. The Eshelby tensor in nonlocal elasticity and in nonlocal micropolar elas- 
ticity. Journal of Mechanics of Materials and Structures, 1, 325-337 (2006). 

[25] Gurtin, M. E., Fried, E. and Anand, L., The Mechanics and Thermodynamics of Continua, Cambridge Uni- 
versity Press, Cambridge, 2010. 

[26] Arnold, V.I. Mathematical Methods of Classical Mechanics, Springer, New York, 1978. 

[27] Zaixing Huang. The damage models based on the representation of nonlocal residual. Mathematics and Me- 
chanics of Solids, DOI: 10. 1 177/10812865 1 1415992 (201 1). 

[28] Fletcher, D. C. Conservation laws in linear elastodynamics. Archive for Rational Mechanics and Analysis, 
60, 329-353 (1976). 



